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Abstract. The first-order Fermi acceleration process at an ultra-relativistic shock wave is expected to create 
a particle spectrum with the unique asymptotic spectral index <7 7 »i ~ 2.2 . Below, we discuss this result and 
differences in its various derivations, which - explicitly or implicitly - always require highly turbulent conditions 
downstream of the shock. In the presence of medium amplitude turbulence the generated particle spectrum can be 
much steeper than the above asymptotic one. We also note problems with application of the pitch angle diffusion 
model for particle transport near the ultra-relativistic shocks. 
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1. Introduction 

Ultra-relativistic shock waves suggested to be sources of 
gamma-ray bursts are also expected by some authors to 
produce ultra-high-energy cosmic ray particles. A process 
of the first-order Fermi acceleration in such shocks was dis- 
. cussed in a series of papers by Bednarz & Ostrowski (1997, 
1998; see also Bednarz 2000a, b), Gallant & Achterberg 
(1999; see also Achterberg et al. 2001), Kirk et al. (2000) 
and Vietri (2002). Below, in section 2, we briefly com- 
pare and discuss different approaches to the considered 
acceleration process, leading to the asymptotic spectral 
index er 7> i w 2.2. We note an important fact that in 
order to derive this result, essentially all these studies 
consider the large amplitude magnetic field perturbations 
near the shock, with the turbulence power concentrated 
in the short wavelength range. The particle energy spec- 
tra generated in shocks propagating in a mildly turbu- 
lent medium, with the limited turbulence downstream of 
the shock, can be much steeper in ultra-relativistic shocks 
than the above 'asymptotic' one. 

In the discussion below we neglect the strictly paral- 
lel shocks, where some of our objections can be invalid. 
However, such shocks are not expected to frequently oc- 
cur in the universe. 
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2. On the first-order Fermi acceleration at 
ultra-relativistic shocks 

The first-order Fermi acceleration process at an 
ultra-relativistic shock wave involves extreme particle 
anisotropy at the shock in the upstream plasma rest frame 
(UPF), and more mild distributions in the shock nor- 
mal rest frame or the downstream plasma rest frame (cf. 
Begelman & Kirk 1990). Let us consider an individual 
cosmic ray particle acceleration starting with a particle 
crossing the shock upstream (cf. a detailed discussion by 
Gallant & Achterberg 1999). Then, in UPF, its momen- 
tum is nearly parallel to the shock normal. When the shock 
Lorentz factor is large (7> 1) the particle moves in front 
of the shock for a time required for a slight deflection of its 
momentum allowing the shock to overtake it and transmit 
to the downstream region. The deflection proceeds due to 
the magnetic field upstream of the shock, consisting of 
the large scale smooth background structure perturbed 
by the MHD fluctuations. This tiny change of particle 
momentum upstream of the shock allows for its transmis- 
sion downstream of the shock, where - due to the Lorentz 
transformation with a large 7 - its momentum direction 
can be changed at a large angle with respect to its origi- 
nal direction before the transmission upstream. Such large 
amplitude angular scatterings can enable a finite fraction 
of particles to follow trajectories leading to the succes- 
sive transmissions upstream of the shock. Repeating of 
the described loops, with each roughly doubling the par- 
ticle energy, leads to formation of the power law particle 
spectrum. Several authors (Bednarz & Ostrowski 1998, 
Gallant & Achterberg 1999, Gallant et al. 1999) discussed 
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this process leading to formation of the spectrum with 
the energy spectral index a rj 2.2 at 7 3> 1 . Essentially 
the same results were obtained within different approaches 
presented by the above authors and by Kirk et al. (2000) 
and Vietri (2002). We do not consider here the papers 
postulating the acceleration process, but not discussing 
details of the proposed mechanism. 

The work of Bednarz & Ostrowski (1997, 1998) was 
based on Monte Carlo simulations of particle transport 
governed by small amplitude pitch angle scattering. Thus, 
depending on the scattering parameter At (a mean time 
between successive scattering acts) and AQ max (the maxi- 
mum angular scattering amplitude), we were able to model 
situations with different mean field configurations and dif- 
ferent amounts of turbulence. One should note that the 
mean field configuration downstream of the shock was de- 
rived here from the mean upstream field using the appro- 
priate jump conditions and trajectories of particles inter- 
acting with the shock discontinuity were derived exactly 
for such fields. A particle trajectory was derived in the re- 
spective local plasma rest frame, with the Lorentz trans- 
formation applied at each particle shock crossing. The ap- 
proach takes into account correlations in the process due 
to the regular part of the magnetic field, but irregularities 
responsible for pitch angle scattering are introduced at 
random. In order to model particle pitch angle diffusion 
upstream of the shock, with nearly a delta-like angular 
distribution 6 ~ 7 (6 - a momentum vector inclination 
to the shock normal) , an extremely small scattering ampli- 
tude should be used, AVL max <C 7 _1 0- Increasing the shock 
Lorentz factor results in decreasing the momentum per- 
turbation required for its transmission downstream and 
leaves a shorter time for this perturbation, t\. In the ap- 
plied pitch angle diffusion approach, the momentum vari- 
ation due to the regular component of the magnetic field 

1/2 

scales like t±, whence the diffusive change scales like t x . 
Thus growing 7 leads to decreasing t\ and the diffusive 
term has to dominate at sufficiently large 7. It is the rea- 
son why in our simulations the orientation of the regular 
magnetic field ceases to play a role in the limit 7 — > 00, re- 



1 The condition AQ max <C 7 -1 lead to the excessive com- 
putation times. Thus, in our simulations we used a relatively 
'large' maximum amplitude AQ, max = §7 , but comparison 
to the results obtained with smaller AQ max revealed only in- 
significant differences in the obtained particle spectra (as mea- 
sured at the escape boundary placed at 4r 3 downstream of the 
shock; r g is a particle gyroradius). One should note that due 
to the above choice of AQ max and because in our numerical 
code the relative particle velocity with respect to the shock 
instead of the velocity in the plasma rest frame was improp- 
erly used for the particle weighting function in the plasma rest 
frame, the angular distributions presented by us (Bednarz & 
Ostrowski 1998) are slightly different in comparison to the re- 
sults of Gallant et al. (1999) and Kirk et al. (2000). However, 
this error leads only to a wrongly presented angular distribu- 
tions and it does not influence particle distributions considered 
in simulations and the derived spectra. 



suiting in the spectral index convergence to its asymptotic 
value. 

However, one should note that with decreasing At 
and AQ max , when the interaction proceeds at the sub- 
resonance (<C r g ) spatial scale, a serious physical problem 
with the applied approach appears. In order to scatter 
particle momentum uniformly within a narrow cone cen- 
tred on the initial momentum, it requires the short wave 
turbulence to be non-linear at the shortest scales. In our 
discussion of the 'effective' magnetic field, B e , in the pitch 
angle diffusion simulations (Bednarz & Ostrowski 1996; cf. 
Appendix below) we evaluate the lower limit of such an 
effective field from the curvature of simulated particle tra- 
jectories as 



£ e = 5o^l+(o.67^p) 2 , (1) 

taking into account both the background uniform field B a 
and the turbulent component evaluated with the use of the 
scattering parameters AQ max and At (in this expression 
At is given in angular units, it stands for cAt/r g (B )). 
Assuming the constant pitch angle diffusion coefficient 
(oc (AVL rnax ) 2 / At) for a series of computations involving 
smaller and smaller AQ. max ~ 7" 1 we had to use At, 
which scales like 7~ 2 . As a result, to be consistent with 
the assumed scattering model, for large shock Lorentz fac- 
tors the effective magnetic field increases to large values 
due to the required growing power being concentrated in 
the short wave turbulence, B e oc jB . Such conditions 
seem to be unrealistic at least upstream of the shock. 

An analogous pitch angle diffusion modelling ap- 
pended the considerations of Gallant et al. (1999; for 
a more detailed description see Achterberg et al. 2001). 
They considered the highly turbulent conditions near the 
shock leading to the particle pitch angle diffusion with 
respect to the shock normal, i.e. the regular part of the 
magnetic field was neglected. These computations gave es- 
sentially the same spectral indices as the asymptotic one 
derived by Bednarz & Ostrowski (1998). Also, in a variant 
of this model with uniform magnetic field upstream of the 
shock and fully chaotic turbulent field downstream, the 
resulting spectral index did not vary substantially. The 
physical content of the discussed model is substantially 
different from the Bednarz & Ostrowski one because it 
neglects the influence of the uniform field (or long wave- 
length perturbations with A > r g ) resulting in magnetic 
field correlations at both sides of the shock. Thus, it pro- 
vides spectra with the asymptotic spectral index at quite 
moderate 7 ~ 10, a feature also present in the Bednarz 
& Ostrowski simulations for parallel shocks. However, if 
the amplitude of the magnetic field turbulence is lim- 
ited, these simulations cannot reproduce spectrum steep- 
ening (or flattening at intermediate Lorentz factors) in the 
presence of oblique magnetic fields (cf. Ostrowski 1993, 
Bednarz & Ostrowski 1998, Begelman & Kirk 1990). Both 
the above models describe essentially the same physical 
situation for shocks propagating in the highly turbulent 
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medium and, of course, in rarely - if ever - occurring par- 
allel relativistic shocks. 

An alternative discussion of the acceleration process 
presented by Gallant & Achterberg (1999) was based on a 
simple turbulence model. In their approach a highly tur- 
bulent magnetic field configuration was assumed upstream 
and downstream of the shock, idealized as cells filled with 
randomly oriented, uniform (within a cell) magnetic fields. 
With such an approach, particles crossing the shock en- 
ter a new cell with a randomly selected magnetic field 
configuration. Thus, there always occur configurations al- 
lowing some particles crossing downstream to reach the 
shock again and again. As a result of successive energy in- 
creases of the same finite fraction of accelerated particles, 
the power law spectrum is formed. In this model there is 
no need for upstream magnetic field perturbations if the 
considered oblique magnetic field configuration can turn 
all upstream particles back to the shock. 

Two quasi-analytic approaches to the considered ac- 
celeration process were presented by Kirk et al. (2000) 
and Vietri (2002). Both provide methods to solve the 
Fokker-Planck equation describing particle advection with 
the general plasma flow and the small amplitude scat- 
tering of particle pitch angle as measured with respect 
to the shock normal. The important work of Kirk et al. 
modified the Kirk & Schneider (1987) series expansion 
approach to treat the delta-like angular distribution up- 
stream of the shock. An analytically more simple Vietri 
approach applies convenient ansatz'es for the anisotropic 
upstream and downstream particle distributions, resem- 
bling the Peacock (1981) approach to acceleration at 'or- 
dinary' relativistic shocks. Both methods confirm the re- 
sults of the earlier numerical modelling. A deficiency of 
the above semi-analytic approaches is its inability to treat 
situations with mildly perturbed magnetic fields, on av- 
erage oblique to the shock normal. If considered valid for 
different configurations of the mean magnetic field, these 
models require the large amplitude short wave turbulence 
to remove signatures in particle trajectories of the uniform 
background field or of the long wave perturbations. Thus 
it provides an alternative description of the same physi- 
cal situation discussed earlier with numerical methods by 
Bednarz & Ostrowski (1998) in the 7 — > 00 limit or their 
parallel shock results and all other authors applying small 
amplitude pitch angle scattering simulations at parallel 
shock waves. 

3. Conclusions 

The discussed approaches to the cosmic ray first-order 
Fermi acceleration at relativistic shocks yield consistent 
estimates of the asymptotic spectral index « (2.2, 2.3). 
However, the result is not as universal as one could in- 
fer from the convergent conclusions of different authors, 
because all presented derivations require (explicitly or im- 
plicitly) large amplitude MHD turbulence near the shock. 
Only the Bednarz & Ostrowski (1998) modelling allows 
one to treat - in a simplified way - conditions with 



medium amplitude perturbations of the magnetic field. 
In such conditions particle spectra are expected to be 
very steep at high shock Lorentz factors. The spectra 
considered by Bednarz & Ostrowski flatten at large 7 
due to an implicit increase of the short wave turbulence 
in their model, approaching closer and closer the paral- 
lel shock configuration considered by the other authors. 
Until now the situation with the medium amplitude tur- 
bulence, SB < B Q , has not been studied in the limit of 
large 7, however, from comparison with the results of 
Begelman & Kirk (1990), Ostrowski (1993) and of Bednarz 
& Ostrowski (1998) for intermediate shock Lorentz fac- 
tors, we expect very steep spectra to be formed in such 
conditions. Thus, if the conditions with limited turbulence 
are met at a large 7 shock, it can be unable to accel- 
erated particles to very high energies in the first-order 
Fermi mechanism. On the other hand the 'low' energy 
electrons radiating from ultra-relativistic shocks could be 
accelerated by the non-first-order processes, analogous to 
the ones discussed by Hoshino et al. (1992) or Pohl et al. 
(2001). 

The main deficiency of the approaches applying the 
pitch angle diffusion equation, in particular of our own 
attempt to discuss cases with oblique background mag- 
netic fields, is their limitation to particular, highly tur- 
bulent conditions near the shock. This limitation may be 
significant or non-significant, depending on whether such 
conditions exist downstream of the ultra-relativistic shock. 
In the process discussed by Medvedev & Loeb (1999) such 
short-wave non-linear turbulence is created downstream of 
the shock, in the non-resonant wave-vector range for the 
shock accelerated particles. 
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Appendix: Evaluation of the effective magnetic 
field B e in a numerical code applying the discrete 
small-amplitude pitch angle scattering method 

Let us evaluate the magnetic field components responsible 
for regular, A9 r , and turbulent, A8 t , angular deviations 
of the particle momentum during a single particle propa- 
gation time-step At. For the regular deviation due to the 
mean magnetic field B Q 

rAt e 

A0 r ~ "~ Terr - -BoAt , (Al) 
r g (B ) P 

where r g (B) — eB/(pc), p is the relativistic particle mo- 
mentum and e its charge. In the above estimate the unde- 
fined proportionality factor depends on local particle dis- 
tribution anisotropy. For the turbulent component of the 
magnetic field one can evaluate its lower limit analogously 
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to (Al) by assuming the turbulent field component, B t 
to be uniform at the spatial scale cAt: 



A6 t ~ J** = e -B t At 
r g{ B t) P 



(A2) 



For the considered model involving uniform scattering 
within a narrow cone of the opening angle Afl max <C 1, 
the mean scattering angle equals |Af2 maa; . If the magnetic 
field components B Q and B t are oriented randomly with 
respect to each other, then the effective field modifying 
particle trajectory, B e , can be evaluated as B\ = B^ + Bf. 



With the estimate fB t At 



lAQr, 



one obtains B t 



B {^ASl max /[(eB /p)At]} and he can evaluate B e as 



B e — B 



\ 



1 



2 AQ, r . 



eB^ M 
p 



(A3) 



In this rough estimate we give the lower contribution 
from the irregular magnetic field by assuming that tra- 
jectory perturbations are due to structures of the wave- 
length A ~ cAt. In the case of turbulence power concen- 
trated at shorter waves the required wave power is even 
higher to cause the considered scattering. The turbulence 
power with A > cAt would provide correlations of succes- 
sive scatterings, excluded in the considered model. 
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